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Abstract 

The role of higher order coupling of surface vibrations to the relative motion in heavy- 
ion fusion reactions at near-barrier energies is investigated. The coupled channels equa- 
tions are solved to all orders, and also in the linear and the quadratic coupling approxi- 
mations. Taking 64 Ni + 92 > 96 Zr reactions as examples, it is shown that all order couplings 
lead to considerably improved agreement with the experimentally measured fusion cross 
sections and average angular momenta of the compound nucleus for such heavy nearly 
symmetric systems. The importance of higher order coupling is also examined for asym- 
metric systems like 16 + 112 Cd, 144 Sm, for which previous calculations of the fusion 
cross section seemed to indicate that the linear coupling approximation was adequate. It 
is shown that the shape of the barrier distributions and the energy dependence of the 
average angular momentum can change significantly when the higher order couplings are 
included, even for systems where measured fusion cross sections may seem to be well 
reproduced by the linear coupling approximation. 
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I. INTRODUCTION 



The analysis of the fusion process in terms of the barrier distribution has generated 
renewed interest in heavy-ion collisions at energies below and near the Coulomb barrier 
[1-11]. Early studies of sub-barrier fusion reactions compared data and theory in terms 
of the excitation function of the fusion cross section. It is now well established that fusion 
cross sections at sub-barrier energies may be enhanced by several orders of magnitude 
compared with predictions of a one-dimensional potential model, which is due to coupling 
of the relative motion to nuclear intrinsic degrees of freedom[|TIJ. It has been shown 
that under the eigenchannel approximation, these couplings give rise to a distribution 
of potential barriers. Recently, Rowley, Satchler, and Stelson proposed a method to 
extract the barrier distribution directly from measured fusion cross sections 0. Although, 
strictly speaking, this method has a clear physical meaning only in the limit of sudden 



tunneling, i.e. in the limit of a degenerate spectrum of intrinsic motions [13], the same 



analysis was later applied also to the case where the intrinsic motions have finite excitation 



energies [14]. The excitation function of the fusion cross section has to be measured with 



extremely high precision at small energy intervals in order to deduce meaningful barrier 
distributions. Such data are now available for several systems, and have shown that 
the barrier distribution is very sensitive to the nuclear structure of the colliding nuclei. 
The analysis of the barrier distribution clearly shows the effects of couplings to static 
deformations @, |[ ||, ||, vibrational degrees of freedom || ||, transfer channels ||, ||, and 
multi-phonon states^, |J, as well as the effects of projectile excitations |J, in a way much 
more apparent than in the fusion excitation function itself. 

Another quantity which has recently received increasing attention in the study of 
heavy-ion sub-barrier fusion reactions is the angular momentum distribution of the com- 
pound nucleus [15-23]. As in the case of the fusion cross section, the angular momentum 
distribution is also affected by the coupling between the relative motion and intrinsic 
degrees of freedom. Experimental data show that the average angular momentum of the 
compound nucleus formed in heavy-ion fusion reactions at sub-barrier energies is system- 



atically larger than the value expected using the one-dimensional potential model [15 



Furthermore it has been pointed out that the angular momentum distribution is also 



sensitive to the details of the coupling! 19 



The fusion barrier distribution, the angular momentum distributions and the fusion 
excitation function all reflect the physical processes occurring in fusion. In comparing 
measurement and calculation, the excitation function shows most sensitively the energy 
of the lowest barrier, whilst the others show more clearly the effect of couplings over the 
whole range of the barrier distribution. Therefore a simultaneous presentation of the data 
and theory in all three forms (if possible) allows the most complete comparison of data 
and theory. 

Theoretically the standard way to address the effects of the coupling between the rela- 
tive motion and the intrinsic degrees of freedom is to solve the coupled-channels equations, 
including all the relevant channels. Most of the coupled channels calculations performed 
so far use the linear coupling approximation, where the coupling potential is expanded in 
powers of the deformation parameter, keeping only the linear term. Whilst this approach 
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reproduces the experimental data of fusion cross sections for very asymmetric systems, it 
does not explain the data for heavier and nearly symmetric systems [[[5, |17[ |2T, 22, |23| . 



Thus, it is of interest to examine the validity of one of the main approximations in these 
calculations, namely the linear coupling approximation, and see whether the effects of 
non-linear coupling improve the agreement between data and the theoretical calculations 
for such systems. Even in asymmetric systems, the non-linear couplings might be impor- 
tant to reproduce precisely measured data. 

The effects of non-linear coupling can be easily studied if the excitation energy of the in- 



trinsic motion is very small so as to allow one to use the sudden tunneling approximation [|13 



The experimental data of the excitation function of the fusion cross section as well as the 
barrier distribution for the 16 + 154 Sm, 186 W reactions were analyzed in this manner 0, ||. 
The effects of higher order couplings on barrier distributions in the limit of zero excita- 
tion energy has been discussed by Balantekin, Bennett, and Kuyucak in the framework 
of the interacting boson model(IBM) [j24{| . However, for nuclear surface vibrations the 
excitation energies cannot be neglected in most cases, and one has to solve full coupled 
channels equations. Due to the complexity of such calculations, very few studies have 
addressed the effects of the higher order couplings for the vibrational motion. Esbensen 
and Landowne expanded the coupling potential up to the second order with respect to 
the deformation parameter, and have shown that second order coupling leads to a better 
agreement between the theoretical calculations and the experimental fusion cross sections 
for reactions between different nickel isotopes [^5]. The quadratic coupling approximation 
was applied also to the 58,64 Ni + 92,100 Mo reactions []20|. There it was shown that the 
experimental data of both the fusion cross sections and the angular momentum distribu- 
tions are well reproduced by the coupled channels calculations in the quadratic coupling 
approximation. Coupled channels calculations including full order coupling and the fi- 
nite excitation energy of nuclear surface vibrations have recently been performed for the 
58 Ni + 60 Ni reaction^, p6fl . It is seen that higher order couplings are essential in repro- 
ducing the experimental data for this system, and the shape of the barrier distribution 
changes drastically when the effects of the higher order couplings are taken into account. 
Refs.|7J, do not, however, discuss the quality of the quadratic coupling approximation 
and the convergence of the expansion of the coupling potential. 

Although all the above studies and the multi-dimensional tunneling model in Ref. [27 
show the effects of higher order couplings in specific systems, there has not been any 
systematic effort to identify their degree of importance for different systems. Furthermore, 
it is not obvious whether calculations to all orders are necessary or the expansion up 
to the second order is sufficient. In view of the high precision data that have recently 
become available, a critical examination of the effects of these approximations on the cross- 
sections and barrier distribution is necessary before making quantitative comparisons with 
experimental data. 

In this paper we solve the coupled channels equations including the finite excitation 
energies of the vibrational states, and without introducing the expansion with respect to 
deformation parameters. The results of these calculations for fusion cross sections, average 
angular momenta and barrier distributions are compared with those using the linear and 
the quadratic coupling approximations. The paper is organized as follows. In Sec. II 
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the coupled channels calculations which include higher order couplings are formulated. 
Explicit expressions for the matrix elements of higher order terms in both the nuclear and 
the Coulomb couplings are presented. It is seen that inclusion of up to the first second 
order term in the Coulomb coupling is sufficient, but higher order terms are necessary 
for nuclear coupling. In Sec. Ill the coupled channels equations are solved for the 64 Ni 
_l_ 92,96 2 r systems, where the experimental data of both the fusion cross sections and the 
average angular momenta of the compound nucleus are available. The asymmetric systems 
16 + 112 Cd, 144 Sm are also studied in this section and the calculations are compared with 
measured fusion cross sections, barrier distributions, and the average angular momenta 
of the compound nucleus. The summary is given in Sec. IV. 



II. COUPLED CHANNELS EQUATIONS 
AND COUPLING FORM FACTORS 



Consider the problem where the relative motion between colliding nuclei couples to 
a vibrational mode of excitation of the target nucleus. For simplicity excitations of the 
projectile are not considered in this section. It is straightforward to extend the formulae 
to the case where many different vibrational modes are present and where projectile 
excitations also occur. For heavy ions, to a good approximation one can replace the 
angular momentum of the relative motion in each channel by the total angular momentum 



J[28. 2S]. This approximation, often referred to as no-Coriolis approximation, will be used 



throughout this paper. The coupled channels equations then read 



h 2 d 2 J(J+l)h 2 Tr . . Z P Z T e 2 

+ n ' + V N (r) + + nhuj - E c 



r\ 7 O ' r\ O 1 ^ v \ / ' 1 Mil i ^n{r) + V Km(r)i(r) = 0, 

l\i ar z 2/j,r z r J ^7 

(1) 

where r is the radial component of the coordinate of the relative motion, /1 the reduced 
mass, and Vn the nuclear potential in the entrance channel, respectively. E cm is the 
bombarding energy in the center of mass frame and huj is the excitation energy of the 
vibrational phonon. V nm are the coupling form factors, which in the collective model 
consist of Coulomb and nuclear components. These two components are discussed in the 
following sub-sections. 



A. Coulomb coupling form factors 



We first consider the effects of higher order terms of the Coulomb component. In 
Refs. p|, E31 it has been reported that the higher order Coulomb couplings are not impor- 



tant in heavy ion fusion reactions. However, Ref . [^] studied only the excitation function 
of the fusion cross section, and did not discuss the barrier distribution. On the other 
hand, Ref. 1 24]] ignored the finite excitation energy of nuclear intrinsic motions, though it 
discusses the effects on barrier distribution. Here we investigate the effects of higher or- 
der Coulomb couplings on both the excitation function of the fusion cross section and the 
barrier distribution and do not ignore the energy of nuclear intrinsic excitations. Initially, 
we consider the case where the target has only a single phonon excitation. 
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The Coulomb potential between the spherical projectile and the vibrational target is 
given by 



t/ f\ [a f ZpZr^ , , Z P Z T e 2 AnZ P e 2\> + 1 1 

where pr is the charge density of the target nucleus and Qa'o the electric multipole operator 
defined as 

Qx>o = J drZ T ep T (r)r x 'Yy (r). (3) 

Eq. (2) uses the fact that the angular momentum for the relative motion does not change in 
the no-Coriolis approximation, and that the associated spherical harmonics are evaluated 
at the forward angle f = 0, leading to the factor \J 2X ' A ^ 1 ■ If we assume a sharp matter 
distribution for the target nucleus and a phonon excitation of multipolarity A, then the 
electric multipole operator is given by|3Cf 



3e v / d(0 )\ a ' \ , ,, sy(A' + 2)(2A + l) ( A A A' \, ,1 ... 
Qy = -Z T (R^) < Qi\o6\ : \i + (— ) ^= M Q Q Q («a«a)vo ,(4) 



up to second order in the surface coordinate a^, where is the equivalent sharp surface 
radius of the target nucleus. In the collective model of surface oscillations, the surface 
coordinates ot\^ are treated as dynamical variables. They are related to the phonon 
creation and annihilation operators by 

"a/, = «o(4/i + (-y tfl A-^), (5) 

where «o is the amplitude of the zero point motion. It is related to the deformation 
parameter (3\ by «o — /^a/v^A + 1 ||31|| and can be estimated from the measured transition 
probability using 

_ 1 4tt l B(EX) T 

a ° v / 2A + l3Z T (4 0) ) A V e 2 ■ W 

This equation is valid if the amplitude of the vibration is small, and the transition operator 
is linearly proportional to a . 

The Coulomb components of the coupling form factors V nm in Eq.(l) are obtained by 
taking the matrix elements of Vc between n- and m- phonon states. Since we assume 
that there exists only the one phonon state in the vibrational excitation of the target, the 
Coulomb coupling form factors up to second order of a are given by 



4°>(r) = V£\r) = ^Wlgi^ao (7) 



f\r) = 2 V ± * )(A ' + 2) < A0A0]A'0 > 2 a\Z P Z T e^ Y 

11 V ; ^ ' 8tt(2A'+ 1) 1 r> 



If there exist two-phonon multiplets, then the formalism becomes much more complicated 
in the case of non-linear coupling. In the case of the linear coupling, it is known that 
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the no-Coriolis approximation enables us to replace the couplings to all the members of 
the two-phonon multiplets by the coupling to a single state by making an appropriate 
unitary transformation O, pSfl. This leads to a significant reduction of the dimensions 



of the coupled channels problem. This property is lost if one keeps higher order terms 
of the Coulomb coupling since the radial dependence of the coupling form factor for the 
Coulomb part explicitly depends on the multipolarity of the nuclear excitation. 

We now apply Eqs. (7) and (8) to fusion reactions between two 58 Ni nuclei, where 
the importance of second order couplings in the nuclear interaction has been reported J25 . 



We take into account the quadrupole vibrational state at 1.45 MeV, and truncate the 
whole space at the one phonon state. The parameters for the nuclear potential and the 



deformation parameter from Ref.f25| have been used. Since at this stage we want to 
investigate the effects of higher order Coulomb coupling, a linear coupling for the nuclear 
interaction has been used, for ease of calculation. The coupled channels equations are 
solved by imposing the incoming wave boundary condition in the inner region of the fusion 
potential. We found that the second order coupling in the Coulomb interaction causes no 
visible effects to the fusion cross section. It changes the fusion cross section by only about 
0.2 % in the energy region we considered, i.e. from about 10 MeV below the Coulomb 
barrier to about 10 MeV above the Coulomb barrier. Fig.l shows the barrier distribution 
(d 2 (E c , m a)/dE^ m ) as a function of the bombarding energy. As seen from the figure, 
second order Coulomb couplings modify only very marginally the barrier distribution. 
Further calculations showed that the situation does not change when the value of the 
deformation parameter is varied within physically plausible limits. Therefore we hereafter 
use the linear coupling approximation for the Coulomb coupling and investigate the effects 
of the higher order terms only for nuclear coupling. The matrix elements of the Coulomb 
coupling form factor in Eq. (1) are now given by 



3 v v 2 {Rff [3 x 



n,m+l 

+ y/n + 15 n 

,m— 1 

). (10) 



2A + 1 r A+1 J^H 



Note that we have defined the multi-phonon channels by taking the appropriate linear 
combinations of the multi-phonon multiplets. As remarked before, this is possible only 
for the linear coupling approximation in the Coulomb interaction. 

B. Nuclear coupling form factors 

In the collective model, the nuclear interaction is assumed to be a function of the sep- 
aration distance between the vibrating surfaces of the colliding nuclei. It is conventionally 
taken as 

ViN) ^ aX0) = ~ l + exp [(r -R P - 4^ - V^4%o)/a] ' (U) 

Volume conservation introduces a small term which is non-linear with respect to the 
deformation parameter a\o in the denominator of the above Eq.(ll). This is ignored for 
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simplicity in the present study. As in the case of Eq.(2) for the Coulomb coupling, here 
we consider the coupling form factor for the forward angle, which is needed to obtain 
the coupled channels equations in the no-Coriolis approximation. We assume a Woods- 
Saxon form for the nuclear potential. The structure of the resultant formulae in this 
subsection, however, remain unchanged for other forms of the nuclear potential. Denoting 
the eigenvalue of a\o by x, the matrix elements of the nuclear coupling form factor read 

Vnmi r ) = I dxu* n {x)u m {x) j== — — . (12) 

1 + exp[(r -R P - 4 0) - ^M^ x )/a] 

u n (x) is the eigen function of the n -th excited state of the harmonic oscillator. The 
conventional nuclear coupling form factor in the linear coupling approximation is obtained 
by expanding Eq.(ll) with respect to a\o and keeping only the linear term. 

The expectation value of the nuclear potential in the ground state is often replaced 
by the phenomenological potential 

V N (r) = m , (13) 



which is assumed to be known empirically]^]. If we take this prescription, the nuclear 
coupling form factor in Eq.(l) is calculated as 

00 -Vn 



poo 

Vnm( r ) = / dxU* n (x)u m 
J — oo 



X 



1 + exp [( r - Rp _ R p _ J^R*!>x)/d 



/OO \/q 
1 + exp[(r - R P - RP - ^^R { T ] x)/a] 

The last term in this equation is included to make the coupling interaction vanish in the 
entrance channel. Eq. (14) represents the form factor which contains couplings to all 
orders. We use these form factors in the next section in order to discuss the effects of 
higher order coupling to vibrational modes of excitation of the colliding nuclei on heavy 
ion fusion reactions. 

Instead of introducing a phenomenological potential given by Eq.(13) as the bare 
potential in the entrance channel, one could use V^P in Eq. (12) as the nuclear potential 
in the entrance channel. The use of Eq.(13) makes it easier to examine the convergence 
of the effects of higher order terms by comparing the results of the calculations in the 
linear and the quadratic approximations and the full order calculations. Notice that V c 



(TV) 
00 



is identical with the potential given by Eq.(13) in the linear coupling approximation. 

III. RESULTS: EFFECTS OF HIGHER ORDER COUPLINGS 

A. Nearly symmetric systems 

We now present the results of our calculations of fusion cross sections, average angular 
momenta of the compound nucleus, and fusion barrier distributions. We first discuss heavy 
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nearly symmetric systems. The experimental data of the average angular momentum of 
the compound nucleus for several systems are summarized in Fig. 5 of Ref. [fT^| . It 
suggests that the conventional coupled channels calculations do not work for heavy nearly 
symmetric systems. We analyse in particular 64 Ni + 92 ' 96 Zr reactions which are typical 
examples where the conventional coupled channels calculations with the linear coupling 
approximation fail to reproduce the fusion cross sections and average angular momentum 
data||21||. Our aim is to investigate whether the failure is due to the linear coupling 
approximation by performing linear, quadratic and full coupling calculations. 

We take into account the couplings up to two phonon states of the quadrupole surface 
vibrations of 64 Ni and 92 Zr, and of the octupole vibration of 96 Zr. We also take their 
mutual excitations into account. The excitation energies of the single phonon states in 
Ni, 92 Zr, and 96 Zr are 1.34, 0.934, and 1.897 MeV, respectively. We assumed the ra- 
dius parameter associated with the coupling interactions to be 1.2 fm in all cases. The 



we 



deformation parameter of 64 Ni was taken to be /?2 = 0.19 |25|| . Following Refs. |I4|, [32 
used /3 2 = 0.25 for the nuclear coupling associated with the quadrupole vibration of 92 Zr, 
while the deformation parameter in the Coulomb coupling interaction was estimated from 
the measured B(E2)| value to be 0.108. The different value for the nuclear deformation 
parameter from that of the Coulomb coupling parameter was required in order to fit the 
angular distribution of the inelastic scattering of 16 from 92 Zr at 56 MeV |[32|| . The defor- 



mation parameter (3% of 96 Zr was estimated from the recently measured B(E3)| value |33 
to be 0.268. We assumed the same value for the deformation parameter for the nuclear 
coupling as that for the Coulomb deformation parameter for this nucleus. The nuclear 
potentials used in this paper are the same as in Ref. PTJ . These modify the empirical po- 



tentials of Christensen and Winther|34]] by setting the range adjustment parameter AR 
to be fm. 

The excitation function of the fusion cross section for these two systems obtained by 
numerically solving the coupled channels equations are compared with the experimental 
data in Figs. 2 and 3 (upper panels). The experimental data, taken from Ref. consist 
only of the evaporation residue cross sections, and do not include fission following fusion. 
The dotted lines are the results in the one dimensional potential model, i.e. without 
the effects of channel coupling. As is well known, the experimental fusion cross sections 
at subbarrier energies are several orders of magnitude larger than the predictions of this 
model. The dot-dashed lines are the results of the coupled channels calculations when the 
linear coupling approximation is used, which are similar to the results of the simplified 



coupled channels calculations reported in Ref. [21]. They considerably underestimate the 



fusion cross sections at sub-barrier energies for both systems. The situation is slightly im- 
proved when the quadratic coupling approximation is used, i. e. when the nuclear coupling 
potential up to the second order of the deformation parameter ||25||, is included (dashed 
lines). However, there still remain considerable discrepancies between the experimental 
data and the results of the coupled channels calculations. When we include couplings to 
all order, we get the solid lines, which agree very well with the experimental data. Dra- 
matic effects of the higher order couplings on fusion cross sections are observed, especially 
at low energies. 

The lower panels in Figs. 2 and 3 compare the results of our calculations of the average 
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angular momentum of the compound nucleus with the experimental data as a function of 
the bombarding energy. It is defined in terms of the partial fusion cross section o\ as 

<i>=EWE^ (is) 

I I 

The meaning of each line in these figures is the same as in the upper panels. We again 
observe that the experimental data are much better reproduced by taking the effects of 
couplings to all orders into account. We thus conclude that coupling to all orders are 
essential to simultaneously reproduce the fusion cross sections and the average angular 
momentum data for heavy (nearly) symmetric systems. This is in agreement with the 
calculations required to fit the barrier distribution for 58 Ni + 60 Ni reaction |7|]. 

B. Very asymmetric systems 

We next consider the effects of higher order couplings for very asymmetric systems 
where the product of the charges ZpZ^ is relatively small. For such systems, the coupled 
channels calculations in the linear coupling approximation have achieved reasonable suc- 
cess in reproducing fusion excitation functions. However, no study has been performed to 
see whether the effects of higher order couplings on the angular momentum distribution 
of the compound nucleus and on the barrier distributions are small. In this subsection 
we re- analyse the experimental data for the 16 O + 112 Cd reaction, for which both fusion 
cross sections and average angular momentum data are available lf22| , and those for the 
16 + 144 Sm reaction, where the the fusion barrier distribution has been extracted from 
the precisely measured fusion cross sections^]. For simplicity in the calculations we ig- 
nore excitation of the projectile in both reactions. These effects will be discussed in a 
separate paper where it will be shown that the octupole vibration of 16 leads to a 



static renormalization of the fusion barrier |36|, [37 1 



In calculating the fusion cross section for 16 + 112 Cd scattering, we include the double 
quadrupole phonon states and the single octupole phonon state of 112 Cd and their mutual 
excitations. The excitation energies are 0.617 and 2.005 MeV for the one phonon states of 
the quadrupole and the octupole vibrations, respectively. The deformation parameters of 
the quadrupole and the octupole vibrations are estimated to be /?2 = 0.173 and (3% = 0.164, 
respectively |22]]. The radius parameter in the coupling interaction is taken to be 1.2 fm. 



Following Ref.[E2|, we use a Woods-Saxon potential whose depth, range parameter, and 



surface diffuseness are V=58 MeV, r =1.22 fm, and a=0.63 fm, respectively. 

The upper panel of Fig. 4 compares the results of the coupled channels calculations of 



fusion cross sections with the experimental data taken from Ref.|22[. Compared with the 
symmetric systems studied in the previous subsection, the enhancement of the fusion cross 
sections is fairly small. This is partly because the product of the atomic number ZpZt 
in this asymmetric system is smaller than the symmetric systems. If we take the linear 
coupling approximation and estimate the coupling strength F at the barrier position rp 
of the bare Coulomb barrier, one finds 

- * ^ l- (16) 



I r 2 2A + 1 r x B +1 
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The coupling strength is thus proportional to the product ZpZ-r. This product is 384 
for 16 + 112 Cd scattering, whilst it is 1120 for 64 Ni + 92 > 96 Zr reactions. The coupling 
strength in this asymmetric system is therefore several times smaller than in the symmetric 
systems even though the values of the deformation parameters are similar. Another reason 
that the enhancement of the fusion cross sections is small in very asymmetric systems is 
the small reduced mass. In the WKB formula for the barrier penetrability, the mass 
parameter appears in the exponent. Hence the heavier the mass, the more sensitive the 
penetrability to a slight change of the potential. Even though the results in the linear 
coupling approximation (dot-dashed line) show a relatively small enhancement of the 
fusion cross section compared with the no-coupling limit, there is still a significant change 
in going to second order coupling, and then to all order coupling. The situation is similar 
for the average angular momentum. Thus even in such cases with low ZpZ?, if data of 
high precision are available, it seems that the linear coupling approximation is inadequate 
to allow quantitative conclusions to be drawn from a comparison of data and calculations. 

The role of higher order couplings in very asymmetric systems can be more clearly 
seen by investigating the fusion barrier distributions. Therefore, we next consider the 
16 O + 144 Sm reaction, for which the effects on fusion barrier distributions of couplings 
to phonon states were shown experimentally for the first time|5|, ||], [TlJ] . The authors of 
Ref. have shown that the fusion barrier distribution for this system is intimately related 
to the octupole vibration of 144 Sm, and that the quadrupole vibration plays only a minor 
role. Accordingly, we ignore the effects of the couplings to the quadrupole phonon states 
of 144 Sm and include only the single octupole phonon state at 1.81 MeV. The deformation 
parameter f3 3 = 0.205 was used as in Ref. 0. The ion-ion potential was of a Woods-Saxon 
form. The depth, radius parameter and surface diffuseness were 105.1 MeV, 1.1 fm and 
0.75 fm respectively, as given in Ref. EH . 

The upper panel of Fig. 5 shows the fusion excitation function from Ref. || , and the 
calculations. The meaning of each line is the same as in Fig.2. As was the case for 16 + 
112 Cd, we observe that the agreement of the theory and experiment appears to be improved 
only slightly by the inclusion of coupling to all orders. The barrier distribution, however, 
reveals significant changes due to the higher order couplings( see the lower panel of Fig.5). 
Note that there exist two barriers in the present two channel problem. Comparing the 
results of the linear coupling approximation ( the dot-dashed line) with those of the all 
order coupling ( the solid line), one observes that the higher order couplings transfer some 
strength from the lower barrier to the higher barrier, and at the same time lower the 
energies of both barriers. 

This can be viewed in a different way by performing the diagonalization of the coupling 
matrix at each position of the internuclear separation to obtain the effective barriers, as is 
done in the computer code CCMOD (23| . Fig. 6 shows these effective barriers for s-wave 
scattering. The meaning of each line is the same as in Fig.5. We observe that higher 
order couplings decrease the energies of both the lower and the higher barriers, consistent 
with the barrier distributions shown in Fig. 5. The higher order couplings also increase 
the width of both potential barriers (Fig.6), leading to narrower peaks in the barrier 
distribution. This then results in the apparent better separation between the two barriers 
seen in Fig. 5. 
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For these asymmetric reactions, the couplings are weak as a result of a combination 
of the small product of ZpZ? and the relatively small deformation parameters. In such 
cases the first order approximation might have been expected to be valid. Despite this, 
the calculations which include couplings to all orders show significant differences from 
first order calculations. It is clear therefore that high precision measurements should be 
analysed using all order couplings even when coupling is weak. 

IV. SUMMARY 

We have shown that in heavy ion fusion reactions, higher order couplings to nuclear 
surface vibrations play an important role. Such higher order terms in the Coulomb cou- 
pling can be safely neglected. Previous work indicated that standard coupled channels 
calculations are not very successful in describing fusion of heavy symmetric systems. We 
have shown that the data can be described well by coupled channels calculations once 
couplings to all orders are included. We found that for the 64 Ni + 92 > 96 Zr reactions, terms 
beyond those in the quadratic coupling approximation result in further enhancement of 
the fusion cross sections at subbarrier energies. The additional enhancement is as large 
as that due to the inclusion of quadratic coupling. The inclusion of the coupling to all 
orders is crucial to reproduce the experimental fusion cross sections and the average an- 
gular momenta. We performed calculations also for the 16 + 112 Cd, 144 Sm reactions 
as examples of very asymmetric systems where the coupling is weaker. It is found that 
in such cases higher order couplings result in a non- negligible enhancement of the fusion 
cross sections and a significant modification of barrier distributions as well as the average 
angular momenta. 

High precision fusion cross section measurements to deduce the barrier distribution, 
and measurements of angular momentum distributions, are designed to study the impor- 
tant couplings in a reaction. The sensitivity to couplings is greatly enhanced by perform- 
ing experiments involving target projectile combinations with a large value of Z P Z T . It 
has been shown in this paper that higher order coupling significantly affects the barrier 
distribution and average angular momentum even for weak coupling cases like 16 + 
112 Cd, 144 Sm with values of ZpZt ~400. Thus, spurious conclusions regarding the nature 
of couplings could be reached if high quality experimental data, particularly for heavier 
systems, are compared with calculations performed only with first order coupling. The 
stage has now been reached when the standard codes of the coupled channels calculations 
should be revised to include coupling to all orders. 
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Figure Captions 



Fig.l:The barrier distribution for fusion between two 58 Ni nuclei. The one phonon state of 
the quadrupole surface vibration is taken into account. The nuclear interaction is treated 
in the linear coupling approximation. The dotted line corresponds to the case where the 
Coulomb coupling potential is also treated in the linear coupling approximation, while 
the solid line takes into account non-linear terms up to the second order. 
Fig.2:Excitation function of the fusion cross section (upper panel) and the average angular 
momentum of the compound nucleus (lower panel) for the 64 Ni + 92 Zr reaction. The 
experimental data are taken from Ref.f2]|. The two phonon states of the quadrupole 
surface vibration of both the projectile and the target are taken into account in the 
coupled channels calculations. The dotted line is the result in the absence of channel 
coupling. The dot-dashed and the dashed lines are the results when the nuclear potential 
is expanded up to the first and the second order terms in the deformation parameters, 
respectively. The solid line is the results of the coupled channels calculations to all orders, 
obtained without expanding the nuclear potential. 

Fig.3:Same as Fig.2, but for 64 Ni + 96 Zr fusion. The two phonon states of the quadrupole 
surface vibrations of the projectile and those of the octupole surface vibrations of the 
target are taken into account in the coupled channels calculations. The experimental 
data are taken from Ref.[[H]. 

Fig.4:Same as Fig.2, but for 16 + 112 Cd fusion. In the coupled channels calculations, 
the projectile is assumed to be inert. The one and two phonon quadrupole states and the 
one phonon octupole state of the target are taken into account. The experimental data 
are taken from Ref. 



Fig.5:Excitation function of the fusion cross section (upper panel) and the barrier dis- 
tribution (lower panel) for 16 + 144 Sm fusion. In the coupled channels calculations, the 
projectile is assumed to be inert, while the single octupole phonon state of the target is 
taken into account. The meaning of each line is the same as in Fig.2. The experimental 
data are taken from Ref.|J. 

Fig. 6: Effective potential barriers for the s-wave scattering of 16 from 144 Sm obtained 
by diagonalizing the coupling matrix. The meaning of each line is the same as in Fig.5. 
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